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Abstract —Recently, several high-resolution parameter estima¬ 
tion algorithms have been developed to exploit the structure 
of strictly second-order (SO) non-circular (NC) signals. They 
achieve a higher estimation accuracy and can resolve up to twice 
as many signal sources compared to the traditional methods 
for arbitrary signals. In this paper, as a benchmark for these 
NC methods, we derive the closed-form deterministic R- D NC 
Cramer-Rao bound (NC CRB) for the multi-dimensional param¬ 
eter estimation of strictly non-circular (rectilinear) signal sources. 
Assuming a separable centro-symmetric R- D array, we show that 
in some special cases, the deterministic R- D NC CRB reduces to 
the existing deterministic R- D CRB for arbitrary signals. This 
suggests that no gain from strictly non-circular sources (NC 
gain) can be achieved in these cases. For more general scenarios, 
finding an analytical expression of the NC gain for an arbitrary 
number of sources is very challenging. Thus, in this paper, we 
simplify the derived NC CRB and the existing CRB for the 
special case of two closely-spaced strictly non-circular sources 
captured by a uniform linear array (ULA). Subsequently, we 
use these simplified CRB expressions to analytically compute the 
maximum achievable asymptotic NC gain for the considered two 
source case. The resulting expression only depends on the various 
physical parameters and we find the conditions that provide the 
largest NC gain for two sources. Our analysis is supported by 
extensive simulation results. 

Index Terms —Deterministic CRB, Cramer-Rao bound, non¬ 
circular sources, rectilinear, DOA estimation. 


I. Introduction 

T HE problem of estimating the parameters of multi¬ 
dimensional (-R-D) signals with R > 1, such as their 
directions of arrival, directions of departure, frequencies, and 
Doppler shifts, has been an extensive research area with 
widely-spread signal processing applications in radar, sonar, 
channel sounding, and wireless communications. Recently, 
various high-resolution parameter estimation algorithms such 
as NC MUSIC El, NC Root-MUSIC 0, NC Standard 
ESPRIT [[3], and NC Unitary ESPRIT 0, 0 have been 
developed to exploit the structure of signals from strictly 
second-order (SO) non-circular (NC) sources 0. The term 
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strictly SO NC (also called rectilinear) is based on the fact 
that the non-circularity coefficient of these signals is equal to 
one 0 - Examples of digital modulation schemes that use such 
signals include BPSK, PAM, and ASK. The aforementioned 
NC algorithms that exploit the non-circularity property are 
known to achieve a higher estimation accuracy and can resolve 
up to twice as many sources 0 compared to the traditional 
methods for arbitrary signals 0. However, the “NC gain” 
achieved by these algorithms from estimating strictly non¬ 
circular sources has so far only been quantified through sim¬ 
ulations. Hence, analytical expressions are highly desirable to 
study the properties of the NC gain under various conditions. 
As deriving a generic formulation for an arbitrary number of 
sources is very challenging, special cases can be considered to 
provide insights towards devising more general expressions. 
Based on the first-order performance analysis framework in 
0, the scenario of a single strictly non-circular source for NC 
Standard ESPRIT was analyzed in 0- It was found that no NC 
gain can be achieved in this case. A first attempt at analytically 
computing the NC gain of NC Standard ESPRIT for two 
uncorrelated strictly non-circular sources with maximum phase 
separation was taken in 0 - Despite the derivation of a closed- 
form expression for the NC gain, the considered assumptions 
in 0 are rather restrictive and do not provide the desired 
comprehensive insights. 

The performance of high-resolution parameter estimation 
algorithms is often evaluated by comparing them to the deter¬ 
ministic (conditional) and stochastic (unconditional) Cramer- 
Rao bounds (CRBs) derived in EQ) and E2, respectively. 
Whereas the stochastic data assumption requires both the 
signals and the noise to be complex Gaussian-distributed, the 
deterministic model assumes that the signals are arbitrary non- 
random sequences while only the noise follows a complex 
Gaussian distribution. Both CRBs are equally recognized in 
the literature. For the data model used to describe weak- 
sense non-circular sources whose non-circularity coefficient is 
between zero and one, a stochastic NC CRB has been derived 
in 132. The follow-up papers 112 and on consider further 
variations of the underlying stochastic model assumption. The 
stochastic NC CRB in fl2ll was derived by extending the 
original Slepian-Bangs formula for circular complex Gaussian 
distributions lfl5ll to weak-sense non-circular complex Gaus¬ 
sian distributions. However, this bound does not apply in the 
case of strict non-circularity, as in the weak-sense case, the 
real part and the imaginary part of the signal can be treated 
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as independent random variables. This is not true for strictly 
non-circular sources, where the real and imaginary parts are 
linearly dependent. 

In this paper, we derive a closed-form expression of the 
deterministic R -D NC CRB for strictly non-circular signals 
impinging on an arbitrary R -D sensor array. The derivation 
is based on the conventional Slepian-Bangs formula, which is 
still applicable due to the complex Gaussian noise assumption. 
Note that our initial contribution in IT6l only states the 2-D 
result without providing a proof and further analysis. Based 
on the devised R -D NC CRB and assuming the i?-D array 
to be separable and centro-symmetric, we show that in the 
special cases of equal rotation phases and full coherence of 
all strictly non-circular signals as well as for a single strictly 
non-circular source, the deterministic R -D NC CRB reduces 
to the existing deterministic R -D CRB for arbitrary signals 
HO). This suggests that no NC gain from strictly non-circular 
sources can be achieved in these special cases. Note that the 
single source case of the R -D NC CRB has been analyzed 
in |5j for a uniform R-D array that contains a uniform linear 
array (ULA) in each mode. Here, we provide a generalization 
of this case to arbitrarily-formed (non-uniform) separable and 
centro-symmetric R -D arrays. Furthermore, the fact that twice 
as many sources can be resolved from the strictly non-circular 
data model is highlighted. 

In our second contribution, we assume 1-D parameter esti¬ 
mation and simplify the derived deterministic NC CRB and the 
deterministic CRB for the special case of two closely-spaced 
strictly non-circular sources captured by a uniform linear array 
(ULA). These simplified expressions are subsequently used to 
analytically compute the maximum achievable NC gain, which 
only depends on the physical parameters, e.g., the number of 
sensors, the SNR, the correlation, the phase separation, and 
the location of the phase reference of the array. The devised 
expression is based on a truncated Taylor series expansion for 
closely-spaced sources. This is, however, the scenario, where 
high-resolution algorithms are primarily applied. Due to the 
fact that the NC gain expression is very general, the properties 
of the NC gain are studied in terms of the above-mentioned 
physical parameters. For instance, it is shown that the NC gain 
is largest if the sources are uncorrelated, the phase separation 
is maximum, and the phase reference is at the array centroid. 
Under these conditions, the two sources entirely decouple and 
do not influence each other. 

The remainder of this paper is organized as follows: The 
data model is introduced in Section [II] In Section [Till the 
derivation of the deterministic R -D NC CRB is provided while 
its analysis is presented in Section EY1 The asymptotic NC 
gain for two closely-spaced sources is analytically computed 
in Section[V] SectionlVllillustrates and discusses the numerical 
results, and concluding remarks are drawn in Section IVIII 

Notation: We use italic letters for scalars, lower-case bold¬ 
face letters for column vectors, and upper-case bold-face letters 
for matrices. The superscripts T , *, H , and -1 denote the 
transposition, complex conjugation, conjugate transposition, 
and the inversion of a matrix, respectively. The Hadamard 
product of two matrices A and B is represented by AqB, the 
Kronecker product is symbolized by A® B, and the Khatri- 


Rao product (column-wise Kronecker product) is denoted by 
A o B. The operator vec {A} stacks the columns of the 
matrix A G C Mx7V into a column vector of length MN x 1, 
the operator Tr {A} returns the trace of the matrix A, and 
diagja} returns a diagonal matrix with the elements of a 
placed on its diagonal. The matrix II m is the MxM exchange 
matrix with ones on its antidiagonal and zeros elsewhere. 
Also, the vector 1 m denotes the M x 1 vector of ones while 
1 MxM is the MxM matrix of ones. Moreover, Re{-} and 
Im{-} extract the respective real and imaginary parts of a 
complex number or a matrix, | • | represents the absolute 
value of a complex number, and E {•} stands for the statistical 
expectation. 

II. Data Model 

Let N subsequent time instants of the measurement data 
sampled on an arbitrary separable R -D grief] of size M \ x... x 
Mr be represented by a linear superposition of d undamped 
exponentials in additive noise. The f-th time snapshot of the 
observed samples can be modeled as El 

d R ■ (r) 

X m - L ,...,m R {t) = ^2 Si(t) e ikmrll i +n mi ,...,m R (t), (1) 

i= 1 r= 1 

where m r = 1 ,...,M r , t = 1 Si(t) denotes the 

complex amplitude of the i-th undamped exponential at time 
instant t, and k mr defines the sampling gricfl Moreover, /i]' ! 
is the spatial frequency in the r-th mode with i = 1,... ,d and 
r = 1 ,... ,R, and n mi ,... iT71K (f) contains the additive zero- 
mean circularly symmetric complex Gaussian noise samples 
with variance erjj. 

In the array signal processing context, each of the R -D 
exponentials represents a narrow-band planar wavefront from 
stationary far-field sources and the complex amplitudes Sj(t) 
are the source symbols. The objective is to estimate the Rd 
spatial frequencies /Zj = ■ ■ ■, n[ R ^] T , Vi, from ([]}. We 

will also use the notation fi^ = • ■ ■, //,'/ '] T , Vr, for the 

spatial frequencies of all sources in the r-th mode. 

In order to obtain a more compact formulation of ©, 
we collect the observed samples into a measurement matrix 
X G C MxN with M = Jlf=i by stacking the R spatial 
dimensions along the rows and aligning the N time snapshots 
as the columns. Subsequently, X can be modeled as 

X = AS + N G C MxN , (2) 

where S G C dxN is the source symbol matrix and N G 
C MxN contains the noise samples. Furthermore, A = 
[a(fii ),..., a(fj.d)} £ C Mxd is referred to as the array 
steering matrix, which consists of the array steering vectors 
defined by 

a{m) = a (1) (/4 X) ) ® ® a (R) (// R) ) e C Mxl , (3) 

1 An R -D sampling grid is defined to be separable when it is decomposable 
into the outer product of R one-dimensional sampling grids (8). 

2 The number km r represents the coordinates of the sampling grid along 
the r-th mode. In terms of the spatial domain, it represents the sensor positions 
of the array in a A/2 sampling grid. For a uniform sampling grid, we have 
km r — TTL'p 1 . 
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where a^ r \pP) £ C MrXl is the array steering vector of 
the z-th spatial frequency in the ?’-th mode. An alternative 
expression of A is given by 

A = A {1) oA {2) o---oA { r) , ( 4 ) 

where A^ = \a^ r \p{ a^(pP)] £ C MrXd repre¬ 

sents the array steering matrix in the r-th mode. 

The non-circularity of a random variable can be defined 
through the non-circularity coefficient (6j. For every complex 
random variable Z with zero mean, the non-circularity coef¬ 
ficient is given by 

" = fflfn = °- w - 1 - <5) 

The cases |k| =0 and o < |«| < i represent a circularly 
symmetric random variable and a weak-sense non-circular 
variable, respectively. The case |k| = 1 represents a strictly 
non-circular (strict sense or rectilinear) random variable. It can 
be shown that for |k| = 1, the real part and the imaginary part 
of Z are linearly dependent (6), i.e., Ci • Re {Z} = C 2 ■ Im {Z} 
for constants ci, C 2 £ R. 

In the context of array processing, the assumption of strictly 
SO non-circular sources requires that the complex symbol 
amplitudes of each source lie on a rotated line in the complex 
plane. This scenario is found, for instance, when real-valued 
data is transmitted by distinct sources causing different delays 
that result in different phase shifts. In this case, the symbol 
matrix S can be decomposed as 11 

S = *S 0 , (6) 

where So 6 M' ixJV is a real-valued symbol matrix and \F = 
diagcontains stationary complex phase shifts on its 
diagonal that can be different for each source. 

Using ©, the model (0 can be written as 

X = AVS 0 + N. (7) 


III. Deterministic R- D NC Cramer-Rao Bound 

In many applications, estimating the parameters of R -D 
signals with R > 1, can be of high importance. As a 
benchmark of such estimators, the corresponding CRBs for the 
multi-dimensional parameter estimation case are desirable. In 
this section, we first review the R -D CRB for arbitrary multi¬ 
dimensional signals and then derive the i?-D NC CRB for 
multi-dimensional strictly non-circular signals. Additionally, 
we provide simplified expressions of the respective CRBs for 
the 1-D parameter estimation case. 


case was derived in Col. Its closed-form expression is given 
by 

C= ^ .Re{(D H niD) ®R ( s r) } el“x“, (8) 

where 

n^= I M - A(A H A)~ 1 A u <eC M xM (9) 

and 


d=[dW ••• D( fi )]eC MxM (10) 

with D = [d^\ ..., d^p] £ C Mxd , r = 1,..., R, contains 
the partial derivatives of A with respect to the components of 
i = 1..... f/, i n the r-th mode. The vectors d[ r are given 
by d,P = da{jjLi)/dfi'p, Vi. Writing oq instead of a(pi) to 
simplify the notation and using 0, we obtain 


,(r) (1) 

d) = a) 


® rf„ ( ' 




AR) 


( 11 ) 


where d^f' 1 = da'p Moreover, R'p = Irxr 0 
Rg contains the estimated signal covariance matrix Rg = 
**Rso*> where the real-valued sample covariance matrix 
Rg 0 is given by Rg 0 = S q Sq/N. Note that Rg 0 can be 
written in matrix form as 


Rs 0 


R\ " pi,2p A-A 
P2,iVAA A 


pi,aV A A 
P2,dV A A 


_Pd, iV A A pd, 2 A A A A 

where A = ||soi II 2 /N is the empirical source power of the i- 
th source and sj. £ R lxJV is the z-th row of So- Furthermore, 
the empirical correlation coefficients pij that represent the 
empirical correlation between the z-th and the j-th source 
vector are defined by 







( 12 ) 


Note that Rg 0 is symmetric such that pij = p^i. 

In the special case of 1-D parameter estimation, the array 
steering matrix A reduces to A = [a(pi ),..., a(pd)\ G 
C Mxd and the deterministic CRB matrix in 0 simplifies to 

C=^-R e {(D ii niD)QR r s y 1 £R dxd , (13) 
where D becomes 


D = [di ••• d d ] £ C Mxd (14) 

with di = da{pi) / dpi 1 V i. 


A. Deterministic R-D Cramer-Rao Bound 

In the case of arbitrary signals, the set of parameters that 
needs to be considered for the deterministic R-D CRB is 
given by the angular parameters p, = [/rl 1 ! ,..., ] T £ 
R Rd x 1 , the real part and the imaginary part of the symbols 
s = vecjS 1 } £ C Ndxl , and the noise power er^. For this 
parameter set that contains a total of ( 2N+R)d+1 parameters, 
the deterministic CRB matrix in the f?-D parameter estimation 


B. Deterministic R-D NC Cramer-Rao Bound 

In contrast to the case of arbitrary signals, the set of 
parameters for the strictly non-circular source model in 0 is 
given by the angular parameters p £ R Mxl , the real-valued 
symbols so = vecjS^o} £ R Ar<ixl , the rotation phase angles 
ip £ R dxl , and the noise power a 2 . Thus, the number of 
parameters is now equal to (N + R + l)d + 1, which requires 
the derivation of a new CRB for this parameter set. 
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The resulting closed-form expression for the deterministic 
NC CRB matrix C ( W in the R-D case is stated in the 
following theorem: 

Theorem 1. The R-D deterministic NC CRB matrix C©© for 
strictly non-circular sources is given by 

c(nc) = m ' { ( G2 “ GlG o lG ?) 0 4 fl) 


(GiGq 1 H 0 ) © R[ 


(R) 

So 


So 

rT /~i—l ; 


{Go-HjGo'Ho)®#** 


So 


[Hi - HfG^GT) 0 R (Ii) 

-1 


'So 




Go © 4 ? 
Go© 4 f 


h^g^gJ)qr[ r) 


So 


HiQR 


{Hq Gq 1 H 0 ) © R (r) 


So 


{G 0 - HjG q 1 H 0 ) © rW 


So 


{HjG^Gj) © rW 


iTi©4f 

h!qr { s r) 


-1 


Go-HjG^Ho)®^ 

)RdxRd 


-1 


v s 0 


(15) 


/n / r>\ 

where Rg 1 = 1 RxR © Rs a an d the matrices G n and 
H n , n = 0,1, 2, are defined as 

G 0 = Re{^*A a A^} eR dxd , (16) 

Ho = Im{'i'*A H A’®'} £ R dxd , (17) 

Gi = Re {{Ir © M/*).D H A\I/} £ R Rdxd : (18) 

H 1 =Im{{I R ^^*)D u A^} 6l“ xi , (19) 

G 2 = Re { (I R © **)D k D{I r © M/)} € R RdxRd (20) 
and A and D are given by 0 and (110b . respectively. 

Proof: The proof is given in Appendix [A] ■ 

It should be highlighted that the assumption of the R-D 
array to be separable is not required for the derivation of O 
in Appendix [A] In fact. (IT5l i is valid for arbitrarily formed 
R-D arrayfl where the columns of A and D are represented 
accordingly. However, the separability assumption simplifies 
the further analysis and helps with the presentation of our 
results in the following sections. 

In analogy to the 1-D parameter estimation case of the CRB 
for arbitrary signals in (IT3l) . the deterministic 1-D NC CRB 
matrix is stated in the corollary: 


Corollary 1. The deterministic 1-D NC CRB is given by (IT5l) . 
where Rg 0 reduces to Rs 0 and G i, H i, and G 2 simplify to 


Gi = Rc {T'*£> h AT'} £ R dxd , 

(21) 

iTi = Im {’®'*U h AT'} £ R dxd , 

(22) 

G 2 = Rc £ R dxd 

(23) 


with A and D being defined in (113b and in dm 


array is assumed to be separable and centro-symmetric for 
simplicity. Specifically, we investigate the two cases of equal 
rotation phases and full coherence for an arbitrary number 
of strictly non-circular signals before we focus on the single 
source case (d = 1). It is shown that in these special 
cases, the deterministic R-D NC CRB reduces to the R-D 
CRB. Furthermore, we also analyze the maximum number of 
resolvable NC sources. 

For our analysis, we first refine the model in ([7]). Assuming 
the R-D array to be centro-symmetric, i.e., it is symmetric 
with respect to its centroid, its array steering matrices A©) 
from 0 satisfy lfT8l 

n M7 A (r) * = A (r) A 0) Vr, (24) 

where A 0 '* £ (r dxd j s a unitary diagonal matrix that depends 
on the phase reference. If the r-mode array centroid is chosen 
as the phase reference GE we have A 0 '* = I,i- The phase 
reference S' 1 ' 1 along the r-th mode can be defined by 

M r 

{W = Wr © < 25 > 
m r =1 

Note that d©) is a property of the array and independent of 
If ' 1 . Using (l24l) and (1251) . we can decompose the r-mode array 
steering matrix A©) with an arbitrary phase reference along 
the r-th mode as 


A 0 ) = AWAW. (26) 

The matrix A £ C MrXd j s ,-jje arra y steering matrix whose 
phase reference is located at the centroid of the r-th mode such 
that for A ©), from (l24l) . the identity A ©) = IlM r a4©)* holds. 
Furthermore, the diagonal matrix A©) = dingle© 4 )/2i } i _ 1 
contains the shifts of the phase reference for each By 
inserting (l26b into (l24b . we can easily establish the relation 
A 0 = A©) A©) . Thus, if the actual phase reference is 
at the array centroid of the r-th mode, we have <$©) = 0, 
A©) = Ac = Id, and consequently A©) = A©). 

Based on (l26l) . we can rewrite A in 0 as 

A = (A«A©)) o (A©) A©)) o ■ • • o (AW A («)) 

= A A, (27) 

where A = A 1 -© o A©) o • • • o AW g C Mxd and A = 
A©) • A©) • . . . • AW g £dxd_ 

Inserting (l27t into the expression for strictly non-circular 
signals in (0, we obtain 

X = AA'S'S'o + N = A<f>S 0 + IV, (28) 


where we have defined $ = AM/ = diag |e 3 © i+<5 ©} i _ 1 
Si= 


with 


IV. Analysis of the Deterministic R-D NC CRB 

In this section, we discuss interesting special cases and 
properties of the derived R-D NC CRB, where the R-D 

3 These also include non-separable arrays such as cross-arrays and L-shaped 
arrays. 


A. Sources with Equal Phases 

An interesting special case of the model (l27t occurs when 
the phase references in each of the R modes coincide with the 
centroid of the R-D array, i.e., d©- 1 = 0 Vr such that A = Id 
and A = A, and, at the same time, the rotation phase angles 
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for all d sources are the samfl i.e., ipt = ip V i. Hence, we 
have 

$ = V = e^Id. (29) 

Under these assumptions, the matrices G n , n = 0,1,2, can 
be expressed as 

G 0 = Re {e~ iv I d A H AI d (jv} = Re {a h a} = A H A 

G\ = Re {e~'^I Rd D^AI d ^} = Re {G h A} = D H A 
G 2 = Re {e~^I Rd D H DI Rd ^} = Re {D h D} = D H D 

while the matrices H n evaluate to zero. The proof that the ma¬ 
trices A H A G R dxd , D h A G R Rdxd , and D H D G R RdxRd 
are real-valued can be found in Appendix [B] 

Using these observations, all terms in (fTTt containing Hq 
or Hi vanish and the R-D NC CRB matrix simplifies to 

C(nC) = • { ( G2 - GlG ° ^ 0 < } } _1 

= ^-{( dHti a d )&R { s R)T } 

= C, (30) 

where we have used the fact that Rg^ = R^ = R-s^ for 
\F = e 1 * 1 I d . From ( 130b . it is evident that the R- D NC CRB 
reduces to the R-D CRB if the phase reference is at the R- 
D array centroid and the rotation phase angles of the sources 
are equal. This suggests that no gain from strictly non-circular 
sources can be achieved in this case. 


B. Coherent Sources 


In this section, the case of full coherence is discussed, i.e., 
the correlation coefficients /5,;j between all pairs of sources 
are given by \pij\ = 1 V i,j. For simplicity, we assume 
that all the sources have unit power, i.e., Pi = 1 V i. Under 
these assumptions, the sample covariance matrix takes the 

A f D\ 

form R So = 1 dX d such that Ry = l RdxRd . Hence, all the 
Hadamard products with Rg^ in the R-D NC CRB matrix in 
fl3 can be omitted and the remaining parts are arranged in 
the following form 


• c (nc)_1 = G 2 - G± Go 1 + GT 1 FfoG _1 JTo GT 1 


2 N 


+ ifiGn |Gn 1 + G^HoG-'H^Go 1 


G 


+ GiGq 1 H q G~ 1 hJ - H 1 G~ 1 Hj 


= G 2 - (g, - H 1 Go 1 Ht) (g 0 - HqGq 1 Hq^ 

- (H 1 - GiGq 1 H q ) (g 0 - H T 0 G^Ho)~ l Hj, 


G? 

(31) 
Gj 

(32) 


where in OTI l we have defined G = Go — HqGq 1 Ho and 
replaced the terms in the square brackets by applying the 
converse of the matrix inversion lemma, yielding the matrix 


4 The same behavior applies to the more general case of equality modulo 
7T, i.e., (fi = <p + ki • 7r, ki E Z for i = 1,2,... ,d. For simplicity of 
presentation, we assume the angles to be equal, this generalization is however 
straightforward. 


(Go-HoG^Hj)-\ Note that (l32l) can be transformed into 
the block matrix form 


2 

y-Y(nc )^ 1 

2 TV ' 


= G 2 


[Gi 


H i] 


'Go 

Ho 

-1 


H o 

Go. 


kJ 


which represents a very interesting simplification of the orig¬ 
inal expression. 

In the next step, we rewrite the R-D CRB matrix for arbi¬ 
trary signals in ® for the case of coherent sources in a similar 
form. Under the aforementioned assumptions, the R-D sample 

A, / D\ 

covariance matrix is given by R' s 1 = l RxR ®(' i &*l dxd ty) = 
( I R 0 i f’*)l RdxRd (I R 0 \F). Hence, we simplify the original 
form of the CRB in (fl3l > into 


^•G-^Re^niD) 


0 ((I* ® **)l RdxRd (I R ® *))} 

(33) 

= R,e{(I K ® V*)D h D(I r ® \F) — (J fl 0 \F*).D h 

A\F 

• (^*A a A^)~ 1 ^*A H D(/ fl 0 'F)} 

(34) 

= Re{G 2 +jTf 2 -(Gi+jifi) 


• (G 0 +jJTo) _1 (Gi — 

(35) 


where we have introduced additional matrices \F in ( 1341 by 
noting that \F\F* = I d . To proceed we require the following 
lemma: 


Lemma 1. The inverse of a full rank complex-valued matrix 
C = A + j B G C nxn with the real part A G R" xn and the 
imaginary part B G R" xn can be split into its real part and 
its imaginary part as follows: 

G _1 = (A + BA^B)" 1 -j A-'B [A + BA~ 1 B) 1 
if A and (A + BA~ 1 B^j are invertible. 

Proof: To prove this lemma, it is sufficient to multiply G 
with C 1 and show that the result is the identity matrix. ■ 

Applying Lemma □ to 63, we split G 1 into its real 
and imaginary part. After some elementary operations and 
using the fact that = — Ho, we obtain equation (l32l > and 
consequently, we have = G. Thus, both R-D CRBs 

become equal if all the sources are coherent. Note that this 
result is valid for arbitrary R-D arrays as the assumptions 
of separability and centro-symmetry were not used in the 
derivation. Analogously to the special case considered in the 
previous subsection, our findings suggests that no NC gain can 
be achieved for coherent sources. 


C. Single Source Case 

The expression for the deterministic R-D NC CRB is 
formulated in terms of the matrices A, \F, D, and the sample 
covariance matrix Rg 0 . Consequently, it provides no explicit 
insights into the parameters of physical significance, e.g., 
the number of sensors M, the correlation coefficient p, the 
source separation. Knowing how the CRB scales with these 
parameters can facilitate array design decisions on the number 
of required sensors to achieve a certain performance under 
specific conditions. As establishing a generic formulation for 
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an arbitrary number of sources is very challenging, we can, 
however, consider special cases such as the single source case. 
Note that this scenario of the R -D NC CRB has been analyzed 
in 0 for a uniform R -D array containing a uniform linear 
array (ULA) in the ?’-th mode and we found that no NC gain 
can be achieved in this case. Here, we provide a generalization 
of the previous results and simplify the R -D NC CRB for non- 
uniform centro-symmetric and separable R -D arrays. 

So far, we have shown that, from the i?-D NC CRB, no 
NC gain can be obtained if the sources have the same rotation 
phase while the phase reference is at the array centroid or 
if the sources are coherent. As the single source case is a 
special case of each of these two properties, i.e., \F = or 

A f D\ 

W So = lftx r, we can directly conclude that the R -D NC 
CRB and the R -D CRB must be equal for this case as well, 
which is in line with our results in 0. 

The simplified expression of the deterministic R -D NC CRB 
for a single strictly non-circular source is shown in the next 
theorem: 


Theorem 2. For the case of an x ... x Mr (M-element) 
separable R-D array with 5 V) = 0 V r, i.e., the phase 
reference of the centro-symmetric array is at the centroid, and 
a single strictly non-circular source (d = 1), the deterministic 
R-D NC CRB can be simplified to 


G (nc) = diag { [c (nc)(1) ,..., C (nc)(ii) ] } 


G R RxR (36) 


with 


£j(nc)( r ) 


1 M r 1 

Q 2 M Y Mr k 2 r ’ 

^ Z_-/m r =l ^m r 


(37) 


where q represents the effective SNR g = 
being the empirical source power given by 
and Sq G R 7Vxl . 


NP/cr 2 with P 
P = Nil l/N 


Proof: The proof is given in Appendix O ■ 

For the special case of a uniform R-D sampling grid, the 
R-D NC CRB expression from Theorem [2] is simplified in the 
following corollary: 


Corollary 2. For an M-element uniform R-D array with an 
M r -element ULA in the r-th mode and a single strictly non¬ 
circular source (d = 1), the deterministic NC CRB for the r-th 
mode in 03 can be explicitly expressed as 


where k mr 


^l(lic) ( r ) 

-(M r - 


q ' M(M 2 - 1) ' 

l)/2, ■ ■ •, (M r — l)/2. 


(38) 


The expression (j38j is in line with our previous develop¬ 
ments in 0. Moreover, (l38l > is equivalent to the result for the 
single source case of the deterministic R-D CRB for arbitrary 
signals derived in ED- This fact proves our previous claim 
that no improvement in terms of the estimation accuracy can 
be achieved for a single strictly non-circular source. 


D. Maximum Number of Resolvable Sources 

In the case of arbitrary signals, it is well-known from m 
that the upper limi{3 of sources that can be resolved with M 

5 This limit is not reached with all array geometries. An example for an 
array, which can achieve this limit is a ULA. 


sensors is d = M— 1. However, if the sources are strictly non¬ 
circular, we can estimate the DOAs of even more sources than 
sensors available. In this section, we establish the conditions 
under which the deterministic NC CRB is valid for d > M. 

Firstly, it is not difficult to see that the matrices G n and 
H n , n = 0,1,2, can have a rank larger than M. For example, 
the matrix Go can be rewritten as 

G 0 = Rc {\F*A h A\F} = 


2 Mxd 


TRe{f} 

JLJL 

ARe {\F} 

Aim {\F} 


A Im {IF } 


From this equation, it can be seen that unless the phase matrix 
\F is equal to (F = diag{ej Vi }f =1 with Lpi = ip+ki-n, hi G Z, 
i.e., all the rotation phases are equal modulo 7r, Go has a rank 
larger than M if d > M. This result complies with the one 
from Subsection IIV-AI For the matrices Gi, G 2 as well as 
H n , n = 0,1, similar forms are easily found. 

Secondly, regarding the additional dependence of the NC 
CRB on the sample covariance matrix Rs 0 , we have proven 
in Subsection IIV-BI that the NC CRB reduces to the CRB if 
the sources are coherent. This suggests that for non-coherent 
sources, the NC CRB is valid for d > M. 

Consequently, we can infer for a uniform linear array that 
if the sources are non-coherent, i.e., 


\Pi,j\< 1 Vi ^ j in 1,2, —, d, (39) 

and the rotation phase angles are different, i.e., 

\<Pi—ipj\pO mod 7T Vi p j in 1,2, ...,d, (40) 

the Fisher information matrix has full rank and is invertible as 
long as d < 2 (M — 1). To support our claim, we provide the 
numerical evaluation shown in Table Q] which suggests that 
the condition d < 2 (M — 1) represents an upper limit on the 
number of sources that is resolvable. Therefore, up to twice 
as many signal sources can be resolved compared to the case 
of arbitrary signals. 


V. Achievable NC Gain for Two Sources 

After establishing that according to the f?-D NC CRB, no 
NC gain can be attained for a single source, the question to 
be studied is what is the maximum achievable NC gain if at 
least two sources are not fully coherent, their rotation phases 
are different, and the phase reference is arbitrary. Although it 
is well known that exploiting the properties of strictly non¬ 
circular sources can provide significant gains in reducing the 
estimation error, so far, the NC gain could only be quantified 
via simulations. In this section, we analytically compute the 
maximum achievable NC gain associated with strictly non¬ 
circular sources. As finding an analytical expression for an 
arbitrary number of sources is an intricate task, we limit our 
analysis to the case of two closely-spaced strictly non-circular 
sources. The R-D CRB for arbitrary source constellations 
tends to infinity when the source separation approaches zero. 
This is not always true for the R-D NC CRB as under certain 
conditions, a finite value is reached. This observation motivates 
us to derive simplified expressions of the NC CRB and the 
CRB for the two source case, which are subsequently used 
to analytically compute the maximum achievable NC gain. To 
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obtain generic expressions in terms of the physical parameters, 
the derivations are based on the model in (l28l >. 

For simplicity, we limit our analysis to the 1-D parameter 
estimation case and assume a ULA composed of M isotropic 
sensor elements, which is centro-symmetric. The phase refer¬ 
ence is located at an arbitrary position. For this scenario, the 
array steering matrix A in model (l28l > simplifies to 

A=[a(p i) a(p d )\ e C Mxd , (41) 

where the steering vectors a/pf), i = 1..... c/, are defined as 


a(ni) 



(M-l) 

2 


Mi 




(M — 1) 
2 



(42) 


After inserting (TTTb into the expression (l2Tb . it is once more 
apparent that if the phase reference is at the array centroid, we 
have 6 = 0 and consequently A = I d . Moreover, if the phase 
reference is at the first element, we have S = (M — l)/2. 


for the single source case discussed in the previous section. 
Hence, the two sources entirely decouple as if each of them 
was present alone. 

Remark 2: Another special case occurs when the two 
sources approach each other, i.e.. A// approaches zero. In the 
CRB for arbitrary sources this always implies that the CRB 
tends to infinity. This is, however, not always true for the NC 
CRB. The limit can be computed as 


lim Tr 

A/i —^0 


(c (nc) ) = — 
i j i — / 


6 


M/M 2 - 1) 

i ei + h 


sin (A</>) QiQ2 


(45) 


Thus, for p < 1 and Acf> > 0, a finite value is reached. If we 
have p = 0 and A <f> = 7r/2, the limit (l45l > corresponds to (l44l) . 
and for p = 1 and A cf> = 0, the limit tends to infinity as the 
NC CRB matches the CRB. 


A. NC CRB for Two Closely-Spaced Sources 

The result obtained by simplifying the NC CRB for two 
closely-spaced sources can be summarized in the following 
theorem: 

Theorem 3. For the case of an M-element ULA (1-D) and 
two closely-spaced strictly non-circular sources (d = 2), the 
deterministic NC Cramer-Rao bound can be simplified to 
expression (143b below. In (143b . we have defined Ap = \p 2 ~Pi I 
and A<j) = Ap + 6Ap with Ap = \p 2 — 'FtI- Moreover, 
Qi = NPi/a i = 1,2 represents the effective SNR of each 
of the two sources. 

Proof: The proof is given in Appendix [D] ■ 

It is worth highlighting that the analytical expression in 
(l43l > is only an approximate result as the derivation involves 
a Taylor series approximation for small Ap, where the higher 
order terms beyond 0(Ap 4 ) have been neglected. Therefore, 
(l43l > becomes accurate if A p is small. 

Also, note that the behavior of the simplified NC CRB 
in m is symmetric in Ap as the two sources can be 
interchanged. Moreover, as any real-valued data stream can 
be multiplied by the factor —1, which represents a phase shift 
of 7T, it is also 7r-periodic. Combining these two results, only 
the interval Ap £ [0,7r/2] must be considered and the general 
behavior of the NC CRB can be extracted from this interval 
by mirroring and periodification. Consequently, the maximum 
phase separation is given by Ap = 7r/2. 

Based on the result in (143b . simplified expressions for 
several special cases can be deduced, e.g., for two uncorrelated 
(p = 0) or coherent (p = 1) sources as well as for Afi = 0 or 
Ac/r = 7r/2. 

Remark 1: One specific case that is worth highlighting is 
the case p = 0 and A </> = 7r/2, where Ap = 7t/2 and 5 = 0. 
Under these conditions, the NC CRB for two sources in (l43l > 
simplifies to 

Tr (C (nc) \ ---• ^ + (44) 

l J M/M 2 — 1) p!p 2 ’ 1 ’ 

which is independent of Ap. As (l44l) resembles the expression 
for a single source in (l36l >. it is apparent that the individual 
NC CRB for each of the two sources represents the NC CRB 


B. CRB for Two Closely-Spaced Sources 

The corresponding expression of the simplified CRB for two 
closely-spaced sources is stated as follows: 


Theorem 4. For the case of an M-element ULA (1-D) and two 
closely-spaced sources (d = 2), the deterministic Cramer-Rao 
bound can be simplified to expression © below. 


Proof: The proof is given in Appendix [E] ■ 

In analogy to the result for the NC CRB, (l46t becomes exact 
for small Ap and the higher order terms beyond O/Ap 4 ) of 
the Taylor series expansion are negligible. 

Again, more simplified expressions for several special cases 
can be derived from (l46b . e.g., p = 0, p = 1, Af = 0, or 
A cf> = 7t/2. 

Remark 3: A very interesting property of the CRB can be 
shown for p = 1 and Af = 7r/2 with 6 = 0. For these 
parameters, we can reduce the CRB in (l46l > to 


Tr{C} 


1 


360 


Ap 2 M(M — 1)(M — 2)(M + 2)(M + 1) 
Qi + £?2 


Qi 02 


(47) 


which corresponds to the expression of the CRB for p = 0 
and arbitrary Acf>. This implies that a rotation phase separation 
of 7 t/ 2 decorrelates two coherent sources. 

Remark 4: In contrast to the NC CRB, the limit for the CRB 
is given by 


lim Tr{C} = oo V p, V Afi. (48) 

A / j ,—^0 

Therefore, the NC CRB for strictly non-circular sources ex¬ 
hibits substantial benefits compared to the CRB if the sources 
are closely-spaced, incoherent, and have a non-vanishing phase 
discrimination Acf>. 


C. Analytical NC Gain for Two Closely-Spaced Sources 
Based on the simplified expressions for the two-source case 
of the NC CRB in (l43l l and the CRB in (l46l) . we can explicitly 
compute the NC gain for two sources as given in ( l49t below. 
As the derivation of (l49l) is based on ( l43l > and (l46t . it becomes 
accurate for small source separations Ap as well. We can now 















analyze the properties of the NC gain expression for different 
values of p, Ap, and S. 

Remark 5: As already established earlier for an arbitrary 
number of sources, the NC CRB becomes equal to the CRB 
if either p = 1 or if A</> = 0, where A p = 0 and <5 = 0. This 
behavior also reflects in the NC gain computed for two strictly 
non-circular sources as it can easily be verified that for these 
parameter values, the expression (l49t evaluates to 7/ nc ) = 1. 
Hence, no NC gain is obtained in these cases. Note, however, 
that if S ^ 0, i.e., the phase reference is not at the array 
centroid, there may be an NC gain even if A p = 0. 

Remark 6 : By analyzing the NC CRB for two closely- 
spaced sources, we have found that for p = 0 and Ac/) = tt/2 
with c5 = 0, the two sources entirely decouple. Evaluating the 
NC gain expression for these parameters leads to 

~ -1 -^- ( 50 ) 

1 Ap 2 (M-2)(Af + 2) 

Thus, this case represents the largest achievable gain for two 
closely-spaced strictly non-circular sources. It is apparent that 
the NC gain in (l50l > decays in proportion to M -2 but increases 
as Ap decreases. 

Remark 7: The limit of the NC gain for Ap approaching 
zero is given by 

lim ?/ nc ) = oo V p, V Ac/). (51) 

/\/ j ,—^0 

Therefore, the NC gain can theoretically approach infinity if 
the source separation tends to zero. 

D. Two Groups of Equal Phases 

This subsection represents a generalization of the case of 
two uncorrelated strictly non-circular sources to two groups of 
equal phases. Let d mutually uncorrelated sources with unit 
power, i.e., Rs 0 = Id , have the phase angles 

Pi = + h ■ 7T or <pi = p^ +ki - it, 

where ki £ Z, i.e., modulo tt there are only two different 
phase angles: < p^\ and p^ 2 1. Without loss of generality, we can 
reorder the sources such that the d\ sources with phase p^ 
are the sources 1 , 2, ..., d\ and the remaining d — d\ sources 
d\ + 1, d\ + 2,..., d have phase p^ 2 \ Thus, the sources fall 
into two groups, where the NC gain depends on the phase 
separation \p I 2 1 — p^\ of the groups. 


Now, in the special case |c/?I 2 ] — p^\ = 7r/2, i.e., the 
phase separation between the two groups is maximum, it is 
straightforward to see that the matrices Go, Gi, and G 2 are 
block diagonal, i.e., they are zero except for the upper left 
d\ x d\ block matrix and the lower right (d — di) x (d — di) 
block. Combining these matrices and using the fact that the 
correlation coefficients are zero, we can show from the joint 
CRB that the two groups decouple, that is, the first d\ sources 
are completely decoupled from the remaining (d— d\) sources. 
This case can provide a significant gain compared to the CRB 
for arbitrary sources if there are closely-spaced sources that 
belong to different groups. 

VI. Simulation Results 

In this section, we provide simulation results to evaluate 
the behavior of the R- D NC CRB and illustrate our analytical 
results. 

A. Behavior of the Deterministic R-D NC CRB 

In this subsection, we compare the root mean squared error 
(RMSE) of the derived deterministic R-D NC CRB (Det NC 
CRB) to the deterministic R-D CRB (Det CRB) and the 
existing stochastic R-D NC CRB (Sto NC CRB) for weak- 
sense non-circular signals from lfl2ll . Moreover, we include 
the R-D NC Standard ESPRIT (NC SE) and R-D NC Unitary 
ESPRIT (NC UE) algorithms Q as well as their non-NC 
counterparts R-D Standard ESPRIT (SE) and R-D Unitary 
ESPRIT (UE) ESI into the comparison. It is assumed that 
a known number of signals with unit power and real-valued 
symbols drawn from a Gaussian distribution impinge on the 
array. 

Fig. Q] illustrates the RMSE over all sources versus the SNR 
for the centro-symmetric 2-D array (R = 2) in Fig. [2] with 
M = 12, where N = 20 available snapshots of d = 3 sources 
with the spatial frequencies p[ 1%> = 0.25, ptp = 0.25, p ^ = 
0.5, pf 1 = 0.5, p[ 3) = 0.75, and p {3) = 0.75, and a real¬ 
valued pair-wise correlation of p = 0.9. The rotation phases 
contained in 5P are given by p\ =0, pi = 7t/4, and P 3 = 
7t/2. It is apparent from Fig. |T| that the NC SE and NC UE 
algorithms perform close to the derived Det NC CRB and that 
all of these outperform the Sto NC CRB from fl2l . 

In Table [j] we analyze the Det 1-D CRB and the Det 1- 
D NC CRB for a varying number of sources d in case of a 
ULA with M = 4, N = 20, and SNR = 10 dB. The spatial 


Tr |G (nc) | « 50400 • [p 2 Ap 2 M{M - 1 )(M - 2 )(M + 2 ){M + 1) ( Ap 2 (M - 3 )(M + 3) • cos 2 (A</>) + 140 ■ sin 2 (A</))) 
+ (1 - p 2 ) M(M - 1 )(Af + l)(l40 • A p 2 {M - 2 ){M + 2) • cos 2 (A</>) + 8400 • sin 2 (A</>))^ • (43) 

Tr {G} « 50400 ■ (p 2 Ap 2 M(M - 1 )(M - 2 ){M + 2)(M + 1)(a^ 2 (M - 3 ){M + 3) • cos 2 (A</>) + 140 • sin 2 (A(/))) 

+ 140 • (1 - p 2 ) A p 2 M(M - 1 )(M - 2)(M + 2)(M + 1) 


Qi Q2 


(46) 
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Fig. 1. Analytical and empirical RMSEs versus SNR for the 12-element 2-D 
array (R = 2) from Fig. [2] and N = 20, d = 3 correlated sources (p = 0.9) 

at /4 1 * = 0.25, pP = 0.25, ^ = 0.5, ^ = 0.5, = 0.75, 

(3) 

fi 2 1 =0.75 with rotation phases <pi = 0, <p 2 = 7r/4, and <p 3 = 7t/2. 


r^n 


liAj 


r( 2 ) 


• • 
• • 
• • 


• • 
• • 
• • 


r(2) 


Fig. 2. 2-D shift invariance for the depicted centro-symmetric 4x3 sampling 
grid, left: subarrays for the first (horizontal) dimension, right: subarrays for 
the second (vertical) dimension. 


frequencies pi, Vi are distributed equally in the interval [—2,2] 
and the rotation phases cpi, Vi are drawn randomly. It can be 
seen that d max = M — 1 for the CRB and dmax = 2 (M — 1) 
for the NC CRB are the largest numbers of d that lead to an 
invertible Fisher matrix, otherwise, the problem is ill-posed. 
Therefore, twice as many sources can be resolved from the 
strictly non-circular data model. 

B. Analytical Results 

In this subsection, we compare the analytical results “ana” 
in (l43l > and (l46l > to the empirical ones “emp” in (fT3l> and 
Corollary Q] obtained by averaging over 1000 Monte-Carlo 
trials. We have d = 2 sources that impinge on a ULA (1- 
D) with the powers Pi = 0.5 and P 2 = 1.5. The symbols Sq 
are randomly drawn from a real-valued Gaussian distribution. 

In Fig. [3] we display the RMSE of the Det 1-D NC CRB 
and the Det 1-D CRB for d = 2 sources as a function of the 
number of sensors M, where the square root of the analytical 


Fig. 3. Analytical and empirical RMSEs versus the number of sensors M for 
d = 2 correlated sources with N = 10, Ap, = 0.1 rad, p = 0.8, Acp = 7r/3, 
<5 = (M - l)/2, Pi = 1.5, P 2 = 0.5, and crj = 0.032. 


TABLE I 

RMSE FOR A VARYING NUMBER OF SOURCES WITH M = 4 


RMSE 

d = 1 

d = 2 

d = 3 

II 

T3 

d — 5 

d — 6 

fi¬ 

ll 

-9 

CRB 

0.02 

0.13 

0.80 

oo 

oo 

oo 

OO 

NC CRB 

0.02 

0.11 

0.12 

0.14 

0.35 

2.93 

oo 


expressions is taken. The source separation is A// = 0.1 rad 
with p i = 0 and p 2 = 0.1, however, the actual positions 
are irrelevant and have no impact on the performance. The 
remaining parameters are given by N = 10, Acp = tr/3, 5 = 
(M— 1)/2, i.e., the phase reference is located at the first sensor 
element, and er 2 = 0.032. Moreover, the correlation coefficient 
p is set to p = 0.8. It is evident that the analytical results agree 
well with the empirical estimation errors and that both CRBs 
perform similarly for large M. 

Fig. [4] illustrates the asymptotic NC gain in (149b for d = 2 
sources as a function of A/j. The number of sensors is fixed 
to M = 15 and we have p = 0, A cp = 7t/2 as well as S = 0. 
For comparison purposes, we have also included the curves 
for the analytical NC gain of NC SE from J9[ for this specific 
scenario. It can be seen that the NC gain expression becomes 
accurate for small A fi and that it is largest when A// goes 
to zero. Furthermore, the NC gain of NC SE is close to the 
maximum achievable NC gain computed from the NC CRB. 


VII. Conclusion 

In this paper, we have presented a closed-form expression 
of the deterministic R -D NC Cramer-Rao bound for multi¬ 
dimensional strictly non-circular (rectilinear) signals. This 


= * 1 + (l40 • (1 - p 2 )M(M - 1 )(M + 1) • sin 2 (A <j>) (bO - A p\M - 2)(M + 2)))/ 

A p 2 M(M - 1 )(M - 2 )(M + 2 )(M + 1)(A p 2 p 2 {M - 3 )(M + 3) • cos 2 (A 0) + 140 • (1 - p 2 cos 2 (A</>))) 


( 49 ) 
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Fig. 4. Analytical and empirical NC gain versus the source separation A/r 
for d = 2 uncorrelated sources with M = 15, N = 10, A tp = it/2, 5 = 0, 
Pl = 1.5, P 2 = 0.5, and = 0.032. 


bound serves as a benchmark for the recently developed 
algorithms, e.g., R -D NC Standard ESPRIT and R -D NC 
Unitary ESPRIT, that exploit the NC structure of such strictly 
non-circular signals and thus outperform the traditional meth¬ 
ods for arbitrary signals. Based on the resulting R -D NC 
CRB expression and assuming the f?-D array to be separable 
and centro-symmetric, we have shown that in the special 
cases of equal phases and full coherence of the strictly non¬ 
circular signals as well as for a single strictly non-circular 
source, the deterministic R -D NC CRB reduces to the existing 
deterministic R -D CRB for arbitrary signals. This suggests 
that no NC gain can be achieved in these specific cases. 
Furthermore, we have simplified the derived NC CRB and 
the existing CRB for the special case of two closely-spaced 
strictly non-circular signals captured by a uniform linear array 
(ULA). With these simplified CRB expressions, we have then 
analytically computed the maximum achievable asymptotic 
NC gain for this scenario. The resulting expression only 
depends on the various physical parameters, e.g., the number 
of sensors, the signal correlation, etc. Additionally, we have 
analyzed the dependence of the NC gain on these parameters 
to find that the largest NC gain is obtained if the two sources 
are closely-spaced, incoherent, and have a non-vanishing phase 
discrimination. 


Appendix A 
Proof of Theorem |7] 

For convenience, we start our derivation by vectorizing 
the R- D NC data model in £7} by using the property 
vec {AXB} = ( B t 0 A) • vec{X} for arbitrary matrices 
A, 13, and X of appropriate sizes. We obtain 

x = vec{X} = {I N ®Ay)s 0 +n€C MNxl , (52) 

where s 0 = vec {So} = [sj(l),..., sJ(7V)] T G R wdxl 
with so(t), t = 1 being the f-th column of So, 

and n = vec {TV} G ^mtvxi -po su j t t jj e deterministic data 


assumption, the signal vector .s 0 is assumed to be determin¬ 
istic and unknown to the receiver, while the sensor noise n 
is zero-mean circularly symmetric white complex Gaussian 
distributed, i.e., E{nn T } = 0. Hence, the observations x 
satisfy the model 

x~CM(u,H) t (53) 


where v = (In 0 AH/) So and S = o^I MN are the mean 
and the covariance of the array output vector x. 

Let us now define the real-valued vector of unknown 
parameters as 


£=[ M t aj p T (54) 


Here, /x = [juW T ,..., /Lt( fl ) T ] T G R Rdxl is the principal 
parameter vector of interest and So G R Ndxl , p G R dxl , and 
<7^ are the nuisance parameters. As the CRB matrix is usually 
computed by taking the inverse of the Fisher information 
matrix (FIM) J, we first need to calculate J. Due to (l53l >. 
i.e., x is Gaussian distributed, the Slepian-Bangs formulation 
ua of the FIM is still valid for the strictly non-circular data 
model in (l52l >. Hence, the Slepian-Bangs formulation of J for 
the parameter vector £ is given by m 


J P , q = Tr 


dip dd q 


2 -Re 



(55) 


p, Q — 1, , (f? + IV + l)d + 1. 


Note that we are only interested in the CRB for fi, denoted 
as C"( rlt i. Therefore, it is sufficient to compute the upper 
left block of J -1 . In order to find J from (l55l >. the partial 
derivatives of v with respect to the parameters of £ can be 
calculated straightforwardly. We have 

^ = (I N ® (D(I r 0 *))) S ( 0 R) G C MNxRd , (56) 

where D is given in (ITOl) and = [(Jfl0So(l)),..., (I R < 8 ) 
S 0 (N))] T G M. NRdxRd with S 0 (t) = diag |s 0 (f)} G R dxd . 
For the remaining parameters, we get 

|^ = (/^0A^)g C MNxNd , 

os 0 

^=j(/ w ®Tf)S 0 G€ MWxd , ^=0 gR" x1 , 


where S 0 = [S Q (1),..., S 0 (N)] T G R Ndxd . Next, these 
results are combined to obtain 


du 


(I N 0 (D(I r 0 ®))) S { 0 R) (I N 0 A\P), 
j (Jjv 0 A'®') S 0 0] G c MjVx[(ii+JV+1)<i+1 ], (57) 


As for the derivative of S with respect to the only non-zero 
term is 


dcri 


= 1 


MN, 


(58) 


such that 


4J=[0 0 0 i MN }ec MNx ^ R+N+1 ^ d+MN l (59) 
1 L J 
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Inserting (l57l > and (l59l) into (f55l) and bearing in mind that we 
are interested in the //-block of J, only the second term of 
(|55T> is of concern. Therefore, we only consider the non-zero 
block J of J, which is given by 


J = 

where 
G = 


Ju 




Jr 

J S( 

Jc 


¥>,30 


JS( 

Jc 




Re{G H G}, (60) 


(In ® (D(I r © ¥))) S { 0 R) (I n ®AV), 

j(I N ®A'S')S 0 \ £<C MNx{R+N+1)d . (61) 


It is easy to see that J = J T . Consequently, only the block 
matrices on and above the diagonal of J need to be computed. 
For the block matrix we obtain 


In the next step, we need to extract the upper left block of 
J” 1 . To this end, we make use of the following lemma: 


Lemma 2. For matrices A G C pxp , B € C pxq , C G C 7 ' xr , 
D G C qxp , E G C qxq , F G C qxr , G G C rxp , H G C rxq , 
and J G <C rxr the upper left p X p block of the matrix 



A 

B 

C 

K = 

D 

E 

F 


G 

H 

J 


(72) 


is given by 

K i :pi i :p ={A- BE x D - BE~ 1 FSf ; 1 HE~ 1 D 
+ BE-'FSe 1 G + CJ~ 1 HE~ 1 D 
+ CJ 1 HE 1 ■FSe 1 HE- 1 D-CSe 1 G^ , 


J nn - 2 

n 


x>{( 


Ir © So(t))(I R © ¥*)£> H 
•£>(/* ©S'XiflSSoW)} (62) 
= ^ • Re| {(Ir © **)D U D(I R © ¥)) 

iV 'j 

® s o(f))(!it © «o(f)) T f (63) 


t=l 


= — • Re<^ ((I R © ¥*)£> h £>(I* © ¥)) 


© | (64) 


2 TV 


= — • Re<^ ((/« © © ¥)) 


2 TV 
0'S 


0| 

© —S 0 Sq ^ | 

(65) 

,(«) 

'So 

£ j ^RdxRd 

(66) 


where G 2 is defined according to (l20b and we have used 
the fact that diag {a} C diag {b} = C © (ab T ) for arbitrary 
vectors a G C M ,b G C w , and a matrix G G <C MxN . In a 
similar manner, the other blocks of J can be computed. The 
results are given by 

Js 0 ,s 0 = —y ■ I N ® Go € R NdxNd (67) 

o TV" 

J v , v = --G 0 GRso& R dxd (68) 

J M = 4^r T (^®G 1 )GR M><M (69) 

J ao , v = ~-(I N ®Ho)S 0 €R Ndxd (70) 

2N 

<W =-fl So ) GR™, (71) 

where the matrices G n and H n . n = 0,1, 2 are given in (fl6l i- 
(IT9t . Note that we have the symmetries Go = Gq , G 2 = Gj, 
and H 0 = -H$. 


where S E = J ~ HE ^F. 

Proof: The proof of Lemma [2] can easily be constructed 
by applying the inversion formula for a 2 x 2 block-partitioned 
matrix If20l to the 3x3 block matrix in (1721 ) twice. ■ 

Applying Lemma[2]to compute the upper left block of J -1 , 
it is straightforward to obtain the expression in (IT5l) . where we 
have 

Se=™-(g 0 - H^G^Ho) © 

<T n ' ' 

BE~ 1 D = ™ ■ (C.G^GJ) © R ( * ] 

BE X F = ~ ■ (GiGf 1 H 0 ) © R™ 

HE ~ lD = -^r ■ (HoGo 1 ^) ©4? 

HE~ 1 F = ^ ■ (HjGo'Ho) 0 • 

This concludes the proof. □ 


Appendix B 

In this section, we prove that for = 0 V r and subse¬ 
quently A = A, the matrices A K A G R dxd , D U A G M . Rdxd , 
and D h D G ]g RdxRd are real-valued. To this end, we make 
use of the following lemma: 

Lemma 3. For two arbitrary non-singular left TV- real matri¬ 
ces X G C MxN and Y G <C MxN satisfying IIX* = X and 
HY* = Y, respectively, the following identity holds: 

y h x = (nr) H njr = y T nnx* 

= (Y H X)*GR M . (73) 


Therefore, to prove that the aforementioned matrices are 
real-valued, we simply show that the matrices A and D are 
left II-real. It is straightforward to see that due to d <r ' 1 = 0 V r, 
this is the case for A (cf. Equation (l24l >). As for the matrix 
D, we utilize the linearity of the differentiation operator and 
obtain 



an a* 

dp 



(74) 


which also renders D left n-real and concludes the proof. □ 
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Appendix C 
Proof of Theorem [2] 

Evaluating the R-D NC CRB in (IT5l) for the special case 
d = 1, the array steering matrix A reduces to aiji), D = 
[dW,...,<£(*)] e C MxR , V = e^, and R So = s^s 0 /N = 
A where So € K JVxl . Moreover, we choose S^ = 0 V r 
for simplicity. Then, dropping the dependence of a on p and 
utilizing the definitions in (0 and (fill , respectively, we have 


a H a 


d« H 


a 




R R 

a< r >V r > = f[M r = M, 

r —1 r— 1 

{[ aW H aW . d^a^ 

p= i 

p^r 

R /M r 

n« (p)H « (p) - f-j e 

p— 1 \ m r =1 

pjir 

[] 0 W H 0 W • dM H d~M 

P= 1 
p^r 



[] flW"aW 


p= i 

P5^r 



M 

W r 


lvi-r 


m r =l 


rM Vr. 


(75) 


0 Vr, (76) 


(77) 


Using the results in (f75l) - (f77l> . the matrices G n and H n , n = 
0,1, 2, simplify to 


Gn = Af, 


Gr = Ho = Hi = 0, 


G 2 = = 


diag | 




}. 


(78) 

(79) 


where in G 2 , the terms d,( ri ' ,H for r\ ^ r 2 evaluate to 
zero due to m. Inserting these expressions into (IT3T i. the 
remaining part of the R-D NC CRB matrix is given by 


C (nc) _ 


diag{[c (nc)(1) ,...,G (nc)(fl) ]} el 


p RxR 


where 


£f(nc)'- r / _ 


(r) a\ M r 


2 NP M V Mr k 2 

r 2jm r =l V 


V r, 


which is the desired result. 


(80) 

(81) 

(82) 

□ 


Appendix D 
Proof of Theorem [3] 

Based on the model in (l28l > after inserting fiB. we start the 
proof by assuming without loss of generality that the phase 
reference is at the array centroid, i.e., A = Id such that A = 
A and d? = df. Using the results from Appendix iBl we can 
write the real-valued matrices A H A, D u A, and L) u D as 


A h A 


M 

a 


a 

M 


D h A 


0 p 

-P 0 


d h d 


T 7 

7 r ’ 


where we have defined a = a^a 2 = a^a\, ft = d^a 2 
dii , and 7 = d\ [ d , 2 = d}}d t . Then, the matrices Go and 
Ho can be written as 


G 0 = Re {'J>*A h A\P} 
Hq = Im {\f , *A H ASP} 


M 

a ■ cos(A ip) 


a ■ cos(Atp) 

M 


0 a ■ sin(A(/j) 

—a • sin(Ay>) 0 


The matrices Gi, Hi, and G 2 can be expressed in a similar 
manner. In order to obtain an expression of the 1-D NC CRB 
that only depends on the physical parameters, e.g, M, p, A <p, 
etc., we approximate the scalars a, /?, and 7 using a Taylor 
series expansion for small source separations Ap = |// 2 — p,\ j. 
Hence, these approximations become accurate for a small Ap. 
Therefore, for a, we have 


(Af —1) 


(Af — 1) 


a 


= E ~ M + jAp • E 


to 


(M-l) 


(M-l) 


Ap 2 


(M-l) 

2 

E 


TO 2 — 


M-^Ap 2 (M 2 -l)+0(Ap*). 


Note that the terms containing odd powers of m evaluate to 
zero. Similarly, in case of a small A \i, we get for f3 and 7 the 
expressions 


(M-l) 


P = -i- E 


m 


gjraA/i 


m= — 

(M-l) 


(M-l) 


-J • 


to • ( 1 + ]mAp — 


Ap 2 , 


^Ap(M 2 - 1) - ^A M 3 (3M 4 - 10M 2 + 7) + Q{Ap% 

(M-l) 


E rn 2 ■ ej mA ^ 


(M-l) 


—j ■ E to 2 • (1 + ]mAp — to 2 — •••) 


(M-l) 


^(M 2 - 1) - ^A,( 2 (3M 4 - 10 M 2 + 7) + Q{Ap A ). 


Finally, with the sample covariance matrix 

A 


Hs 0 = 


pV A A 
pA A A A 


(83) 


and the help of the Taylor approximation terms above, we 
can evaluate the 1-D NC CRB expression in Corollary Q] for 
two closely-spaced strictly non-circular sources. Due to the 
cancellation of relevant terms when using only approximation 
terms of lower order, we also need to consider higher-order 
Taylor approximation term0 for a, /3, and 7 . After some 


6 Here, we used Taylor approximation terms up to the 6th order. 
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tedious calculations, we obtain 


where 




A+A 

AA 


x 0 + xi A fj, 2 + x 2 A// 4 + • • ■ 

2 /i A/i 2 + 2/2 A,u 4 + 2/3 A/r 6 H-' 


(84) 


(85) 


It is apparent that the first term in the numerator and the first 
two terms in the denominator of (l84l) are dominant. Neglecting 
the non-relevant higher-order terms in the numerator and de¬ 
nominator of (l84l > and applying some algebraic manipulations, 
an expression in the form of (l43l > can be deduced. Finally, 
to make the result more general, we consider an arbitrary 
phase reference and substitute A tp by Ac/) to obtain (l43l i. This 
concludes the proof. □ 


Appendix E 
Proof of Theorem |4] 

The proof of Theorem [4] follows the same steps as the proof 
in Appendix [D] Under the same assumptions, we compute the 
matrices A H A, D u A, and D u D in the same way. The dif¬ 
ference is, however, that we evaluate the 1-D CRB expression 
given in ®. Using the same Taylor series approximations 
as before, we obtain a similar expression as (l84l >. Finally, 
neglecting the non-dominant terms in the numerator and the 
denominator, and substituting A</> for Ac p, we arrive at the 
expression in (l46l > to prove this theorem. □ 
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